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C\ ! Abstract 

^ ■ While near-field scanning optical microscopy (NSOM) can provide optical 

images with resolution much better than the diffraction limit, analysis and 

o ' 

O ■ interpretation of these images is often difficult. We present a theory of imag- 



ing with transmission NSOM that includes the effects of tip field, tip/sample 
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coupling, light propagation through the sample and light collection. We ap- 
ply this theory to analyze experimental NSOM images of a nanochannel glass 
(NCG) array obtained in transmission mode. The NCG is a triangular array 
of dielectric rods in a dielectric glass matrix with a two-dimensional pho- 
tonic band structure. We determine the modes for the NCG photonic crystal 
and simulate the observed data. The calculations show large contrast at low 
numerical aperture (NA) of the collection optics and detailed structure at 
high NA consistent with the observed images. We present calculations as a 
function of NA to identify how the NCG photonic modes contribute to and 
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determine the spatial structure in these images. Calculations are presented 
as a function of tip/sample position, sample index contrast and geometry, 
and aperture size to identify the factors that determine image formation with 
transmission NSOM in this experiment. 

PACS numbers: 07.79.Fc, 61.16.Ch, 42.70.Qs, 42.25.Fx, 42.25.Bs, 42.30.Va 
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I. INTRODUCTION 



Near-field scanning optical microscopy (NSOM) is an exciting new class of optical mi- 
croscopies which can provide optical resolution much better than the diffraction limit . 
One realization of NSOM uses an aperture that is much smaller than the wavelength, A, of 
the light as a nearly point-like (on the scale of A) light source. Typically an optical fiber 
is pulled to a 20-100 nm tip and metal coated, leaving a small hole in the metal coating 
at the tip to provide a nanometer-scale aperture. This aperture is then placed very close 
to the sample surface so that light emitted from the aperture does not diffract significantly 
before reaching the sample and superresolution, well below the diffraction limit A/2, can be 
achieved. Components of the light that are strongly localized laterally by the aperture are 
evanescent, decaying rapidly as they move away from the aperture. Tip/sample distance 
can be adjusted to control the contribution from these evanescent tip fields. 

NSOM images can now be readily obtained. A key step to the further development and 
application of NSOM is learning how to interpret, understand, and analyze these images. In 
NSOM, the excitation, the detection, or both can occur in the near field. Strong coupling 
between the sample and the light source/detector, which is not present in far-field optical 
microscopy, will occur in NSOM. For example, in transmission (reflection) NSOM, light 
from the metallized fiber tip couples to the sample in the near field of the tip, while the 
light transmitted (reflected) from the sample is collected in the far field. The NSOM images 
are influenced both by the strong sample/source coupling and by the far-field optics of the 
collection process. To understand these images one must know how the localized source 
field is influenced by the presence of the sample in the near field, how light scatters from 
the sample, and how light is collected in the far field. In this paper, we present a theory for 
imaging with transmission near-field optical microscopy. Our goal was to develop a theory 
capable of describing all parts of the imaging process. A complete theory is necessary to 
clearly identify and separate the contributions made by each step in the imaging process. 

We apply this theory to understand transmission NSOM images ]7| made recently of 
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a nanochannel-glass (NCG) array ||. These images probe the optical mode structure of 
this two-dimensional (2D) photonic crystal ||. The NCG array studied is a 2D triangular 
array of glass rods in a matrix made from a lower-index glass. By heating and pulling this 
structure the lattice constants of the array can be controlled. High quality NCGs with lattice 
spacings on the order of or much less than A can be made. The surfaces that are scanned 
have nearly flat topography. For these reasons, NCGs provide excellent samples for testing 
NSOM. A comparison between theory and experiment will identify which features of NSOM 
are important in this experiment, what determines the observed contrast and resolution in 
the transmission NSOM images, and which features in the images depend on the details of 
the tip field. 

The index variation of the NCG array is periodic in the 2D plane perpendicular to the 
glass rods, so the array possesses the optical modes of a 2D photonic crystal. Because the 
glass rods have the higher index, the distribution of these photonic modes in this 2D plane 
tends to be larger in the glass rods. However the index mismatch between core and matrix 
glasses is small for the sample studied, so significant coupling of these photonic modes among 
neighboring cores occurs as well. A comparison between theory and experiment will show 
what information about the 2D spatial distribution of the NCG photonic crystal modes can 
be learned from imaging with transmission NSOM. 

We briefly describe in Sec. II the experiment |7j to be modeled and the key findings that 
we wish to understand. The theory used to model images obtained by transmission NSOM 
is presented in Sec. III. Each step in the process is modeled, including the tip fields, the 
optical modes of the sample, field propagation from the source through the sample to the 
detector, and collection of the transmitted light by the detector. The model used for the 
nanochannel glass array is also described in this section. In Sec. IV the calculations done to 
simulate the experiments are presented. The results show that the theory provides a good 
description of the observed images. The results allow us to determine which features of the 
near-field optics most affect these experiments. In this section we also present simulations 
for transmission NSOM of samples with much higher index contrast than that of the samples 
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studied experimentally. Stronger coupling of the near field to the sample photonic modes 
is achieved by increasing the index contrast in the NCG. For this case we clearly identify 
those features in the transmission NSOM that arise from strong coupling of the source to 
the sample in the near field and determine how to enhance this coupling. Conclusions are 
presented in Sec. V. 



II. EXPERIMENTAL RESULTS 

The transmission NSOM experiment we model is shown schematically in Fig. [l|. A metal 
coated fiber NSOM tip was placed about 10 nm from the sample surface. Light transmitted 
through the sample was collected in the far field by an objective with numerical aperture 



NA |12] . The nanochannel glass array sample was scanned in x and y at constant separation 
above the tip to produce a 2D image of the NCG. The images were taken with two different 
wavelengths of light (A = 670 nm and A = 488 nm). Polarization of the light leaving the 
fiber tip was controlled by use of fiber paddles. Light was collected for three different NA 
(NA = 0.28,0.55,0.7). 

The nanochannel glass sample studied in the experiment @ was a 2D triangular lattice of 
one glass (channel glass) embedded in another glass (matrix glass) with similar but slightly 
lower index of refraction. The channel glass was cylindrical, approximately 745 nm in 
diameter with center-to-center nearest neighbor separation of 1.07 ± 0.05 fim. The index of 
refraction of the matrix glass was 1.65-1.68 and that of the channel glass was 0.2% to 1.2% 
higher in the visible range. 

The detailed results for this experiment, including the gray-scale images made by trans- 
mission NSOM and line scans from these gray-scale images, are presented in Ref. 0. We 
summarize here the key findings. For low NA of the collection optics (NA = 0.28), the 
gray-scale images exhibit a triangular array of bright, circular spots with approximately 
55% optical contrast between the bright regions, with maximum intensity h, and the dark 
regions, with minimum intensity Id- The optical contrast is the change in intensity normal- 
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ized by the average intensity 2(1^ — Id) / (h + Id)- The center of a spot corresponds to the tip 
below the center of a channel-glass core. Similar results are seen for both A. For intermediate 
NA (NA = 0.55), the images for the two A are different. For the longer wavelength (A = 670 
nm), the bright spots flatten, the light to dark contrast is reduced to 38%, and the overall 
signal level increases. For the shorter wavelength (A = 488 nm), the bright spots becomes 
rings, brightest when the tip is near but not at the channel glass edge, weaker when near the 
center of the channel glass core, and weakest outside the channel glass. The light to dark 
contrast is reduced to 25%. At the highest NA (NA = 0.70) this ring structure persists for 
the shorter wavelength with the contrast reduced to 15%. For the longer wavelength this 
ring structure begins to appear and the the contrast ratio is reduced to 28%. These results 
are independent of the polarization of the tip field and weakly dependent on the aperture 
size and the tip/sample separation in the near field. Typical line scans are shown in Fig. |2| 
for NA = 0.70 with the tip about 10 nm from the sample. The scans were done along x, 
a direction which goes between two adjacent glass cores, as shown in Fig. |I|. Diminished 
contrast and a weak contrast reversal are observed when the tip and sample are widely 
separated by approximately 1 /xm. 



III. THEORY 

To develop a complete theory for imaging in transmission, we must have models for 
the incident tip field and for the electromagnetic modes of the sample. We must correctly 
account for the coupling of the incident field to the sample modes and for the coupling of 
the excited sample modes to reflected and transmitted fields. Finally, we must model the 
collection of the transmitted light. These steps are modeled in the following four sections. 



A. Tip field 



Rigorous, three-dimensional calculations of the fields emitted by a tip with a nanometer 
aperture are difficult and numerically intensive |]6|Jl~3[|. The Bethe-Bouwkamp (BB) model 
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|10|JTT| is a simple approach that is commonly and successfully used |T|,|14[] instead of the 
numerically intensive approaches to model the tip field. In this paper, we use the BB model. 
As our results will show, the BB model provides an adequate representation for the tip-fields 
in the experiments we describe so we did not test any of the models (see for example |T5HT7| ) 
that improve upon the BB model. 

In the Bethe-Bouwkamp model, the near field of a tapered, metal-coated tip is modeled 
by the near field of the light transmitted by a circular aperture with radius a, in a perfectly 
conducting, thin, metal screen when a linearly polarized, plane-wave field is incident normal 
to the screen (see Fig. |l|). Bethe and Bouwkamp calculated the fields, accurate to order 
k Q a, k = 2tc/X, for this problem 50 years ago fllOljrif . Most NSOM experiments are done in 
the range 0.1 < k a < 1.5, where the BB fields give reasonable results. As pointed out by 
Grober et al, the higher-order corrections to the BB result change only the ratio of electric 



to magnetic field amplitudes without strongly affecting the spatial variation of the field |17 

Bethe originally solved the problem in terms of a fictitious magnetic surface-charge den- 
sity, p m , and a magnetic surface- current density, J m , that exist inside the hole in the screen. 
Bouwkamp corrected the J m determined by Bethe. Magnetic vector and scalar potentials, 
A m and m were determined from J m and p m . The electric and magnetic fields transmitted 
by the aperture, E ap and B ap were obtained by taking appropriate derivatives of the vector 
and scalar potentials. Summarizing, 

E ap = V x A m , (1) 
B ap = -ik A m - V0 m , (2) 

where 

A m = - J J m — dA, (3) 



Pm— dA. (4) 



The integrals are done over the area A of the aperture, and r is the distance from the field 

point to the source point in the aperture. For a plane-wave electric field, linearly polarized 

along x with amplitude Eq and incident normal to the screen, the magnetic charge and 

current densities at a source point (x, y) inside the aperture are, 

2yE 
7Y z y/a z — x z — y z 

and 

2ik xyE 

J ™ ~ My/ai-x*-y>' [ ' 

_ 2ik E (2a 2 -x 2 -2y 2 ) 

and 

Jmz = 0. (8) 

The fields transmitted when the incident field is linearly polarized along y can be obtained 
from Eqs. [!]-[§] by applying the appropriate rotations. Transmitted fields for arbitrary inci- 
dent polarization can be obtained from the appropriate linear combination of these two sets 
of fields. 

All of the experimental images were made with the tip axis, taken here to be along the z 
direction as shown in Fig. [I], parallel to the glass rods in the NCG and normal to the NCG 
surface. To model images, we must do calculations for different positions of the aperture 
in an x-y plane at a fixed distance from the NCG and for different distances between the 
tip and the NCG surface. These calculations are most easily done in x-y Fourier space 
because the NCG is periodic in x and y. To do these calculations we need the 2D Fourier 
transform of the tip fields. We first obtain the 2D Fourier transform of the tip fields in an 
x-y plane a distance z from the plane of the aperture, E ap (k x ,k y , z ) and H ap (k x , k y , z ), 
by numerical fast Fourier transform. From the free-space propagation, we can obtain the 
Fourier transform of the tip field at any other x-y plane in free space at a distance z from 
the aperture. For the electric field: 

8 



E ap (k x , k y , z) = exp(ik z (z - z ))E ap (k x , k y , z ) (9) 

where the wavevector k z for propagation along z is obtained from the free-space dispersion 
of a mode with transverse wavevector (k x ,k y ), k z = \Jk 2 — k x — ky. The phase of k z is 
chosen so that the propagating modes ik z real) are moving forward away from the aperture 
and the evanescent modes ik z imaginary) are damped going away from the aperture. A 
similar equation applies for the magnetic field. Translations of the tip field in an x-y plane 
are achieved by translating the phase of the Fourier transform of the field. Shifting the 
tip from (0, 0, z ap ) to (x ap , y ap , z ap ) changes the Fourier components from E^A^, k y , z ap ) to 
exp( i(k x x ap -\- k y y ap )^jJ*j ap (k x , ky, z ap j. 

B. Optical modes of the Nanochannel Glass Array 

Typically, when the modes of a 2D photonic crystal are determined, the band structure 
for the mode frequency, u versus k x and k y , is found for k z = 0. To model transmission 
through an NCG, we need a different set of the optical modes; we need all propagating 
and evanescent modes, all k z , for a given k x , k y , and u. For that reason, we solve an 
eigenvalue problem for k z rather than an eigenvalue problem for u. We begin from the 
Maxwell equations for electric and magnetic fields, E and B, with frequency uo, 

V-B = 

V • n 2 E = 

V x B = -n 2 i(-)E 

c 

VxE = ^)B, (10) 

where n is the local index of refraction in the NCG. Defining F = nE, then 

,\ , Vn 1 
V x E = V x (-F) = -xF+-VxF, 

n n 2 n 

V-n 2 E = V-nF = nV-F + F- Vn. (11) 

From this we have the following Maxwell equations, 
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V-B = 

Vn 

V-F + F = 

n 

V x B = -m(-)F 

c 

V x F- — x F = m(-)B. (12) 

n c 

Using the definition g = Vn/n, the Maxwell equations become, 

V-B = 
V-F + F-g = 

V x B = -m(-)F 

c 

V x F - g x F = ni(-)B. (13) 

c 

To find the modes of the NCG, we treat the NCG as translationally invariant along z. 
The NCG is periodic in x and y, so n and g have the form 



n(x, y,z) = Y^ n ( G ) exp(i(G x x + G y y)) 

G 

g(x, y,z) = J2 g(G) exp(i(G x x + G y y)). (14) 

G 

The field modes must have the form, 

B m «(x, y,z) = J2 B mq (G) exp(t((q x + G x )x + (q y + G y )y + k^z)) 

G 

F mq (a:, y, z) = ^ F mq (G) exp( ? (( g:c + G x )x + (q y + G y )y + fc^)), (15) 

G 

where the G are the 2D reciprocal lattice vectors for the NCG triangular lattice, q is in the 
first Brillouin zone and /c™ q is the eigenvalue for the m-th mode for a given q. 

Define k x = q x + G x , k y = q y + G y and the convolution (n*F)(G) = J2c n(G — G')F(G') 
where the sum is over the 2D reciprocal lattice vectors G' and similar definitions apply for 
other pairs of n or g convolved with B or F. For a given q, the Maxwell equations couple 
field Fourier components with different G. The equations are (suppressing the Fourier and 
eigenmode indices), 
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ik x B y - ik y B x + i( — )n * F z = (16) 

ik x ¥ y - ik y F x — g x *F y + g y *F X — i(^-)n * B z = (17) 

ik x B x + ik y B y + ik z B z = (18) 

ik x F x + ikyFy + ik z F z + g x * F x + g y * F y = (19) 

ikyB z — ikzBy — — i(— )n * F x (20) 

ik z H x — ik x ~B z = — i(— )n * F y (21) 

ik y F z - ik z F y - g y * F z = i(^-)n * B x (22) 

ik z F x -ik x F z + g x *F z = i(-)n*B y . (23) 



The first two equations, Eqs. (16) and (17), come from equations that have no derivatives of 
the z field components and have no dependence on k z . These equations provide homogeneous 
constraints on the solutions to the eigenvalue problem defined by Eqs. (18)-(23). Written 



in explicit eigenvalue form, Eqs. (18)-(23) become, 

k x B z -(^)n*F y = k z B x (24) 

k y B z + (^-)n*F x = k z B y (25) 

k x B x kyBy = k z B z (26) 

(^)n*B y + k x F z + ig x *F z = k z F x (27) 

—(—)n *B X + k y F z + ig y * F 2 = k z F y (28) 

c 

-k x F x + %g x *F X - k y Fy + ig y *F y = k z F z . (29) 



Eqs. (24)-(29) and the two constraints, Eqs. (16) and (17), guarantee that for a given q 
we will have four modes per G. These four modes correspond to two polarizations for each 
k z and -k z . For k z -> -k z , we have (F x ,F y ,F z ) -> (— F x , —F y , F z ) and (B x ,B y ,B z ) -> 
(B x ,B y , — B 2 ), as can be checked by inspection of Eqs. (16), (17), and (24)-(29). 



11 



C. Field transmission 



We determine the transmission of the incident tip field through the NCG by solving the 
standard boundary matching problem for transmission through a film. We assume that the 
total field in the free-space region containing the tip is the incident tip field plus the field 
reflected by the NCG. Explicitly, the reflected field is 

E r (x, y,z) =J2 E r, P ( k ) exp(i(k x x + k y y + k r z z)), (30) 

k,p 

where the sum is over lateral wavevectors k = (k x , k y ) and the two polarizations p. The z 
wavevector k r z is obtained from the free-space dispersion k r z = ^k^ — k 2 x — k 2 y . The phase of 
k r z is chosen so that propagating modes are forward propagating and evanescent modes are 
damped as the mode moves from the NCG sample back toward the tip. The two polarization 
vectors E rjP (k) are chosen to ensure that V • E = for each free-space mode. Reflections 
off the tip are not included. This is the only approximation, other than the use of the 
Bethe-Bouwkamp model for the source field, that we make for the local fields near the tip. 
As our results will show, including multiple reflections is not important for the experiments 
we model. B r is defined by a similar equation. The field in the free-space region containing 
the collection optics is the transmitted field, 

E t (x, y } z)=J2 E *, P ( k ) exp(i(k x x + k y y + k\z)), (31) 

where k\ is determined from the free-space dispersion. The phase of k\ is chosen so that 
propagating modes are forward propagating and evanescent modes are damped as the mode 
moves from the NCG sample toward the collector. The field in the NCG is 

E ncg (x, y,z) = J2 « mq E mq (x, y, z), (32) 

m,q 

where the ct mq are the amplitudes for the modes excited in the NCG. At each NCG surface, 
we have the boundary conditions that the two tangential electric field components, E x and 
E y , and the two tangential magnetic field components, B x and B y , be continuous. These four 
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conditions define the set of matrix equations that are solved to determine the amplitudes 
for the Fourier components of the reflected and transmitted fields and for the excited NCG 
modes. 

When calculating transmission through a film, one must worry about the effects of mul- 
tiple reflections between the sample surfaces that can produce resonant transmission and 
Fabry Perot oscillations. We calculate transmission for thick NCG samples, as used in the 
experiments. These samples are much thicker than A. We also do calculations for films much 
thinner than A. We study these two limits to investigate the effect of sample thickness on 
image formation. We have no problem doing the calculations when the film is thinner than 
A, because transmission resonances do not occur. For the thick samples, we have to be more 
careful. Resonances can be sharp when the sample has a uniform thickness that is much 
greater than A. In real samples these resonances are broadened because the sample thickness 
is not uniform. However, our calculations are done for samples with uniform thickness. Spu- 
rious effects from sharp resonances in the transmission can occur because we include only a 
finite set of reciprocal lattice vectors in our representations of the index and the fields (Eqs. 
H| and |TB|). Typically, we sum over all reciprocal lattice vectors with lengths less than 3-4 



times the length of the fundamental reciprocal lattice vectors. Numerical diagonalization 
to find the eigenmodes becomes very time consuming if we include more reciprocal lattice 
vectors. Small changes in the cutoffs make only small changes in the modes. This can result 
in large changes of the specific modes that are resonant at a particular sample thickness. 
We have checked that the effects of these resonances are spurious by varying cutoffs for the 
wavevectors we include and by varying the sample thickness. We can eliminate these spu- 
rious effects by increasing the cutoffs. However, for typical cases, it is computationally too 
expensive to increase the cutoff enough to eliminate all spurious effects. To eliminate these 
spurious effects simply, we add a damping factor to all of the propagating NCG modes. For 
all NCG modes with real k z , we include an imaginary term (typically 0.01 (/im) _1 ) to the 
z wavevector to damp waves as they propagate across the sample. This damping reduces 
the effects of the multiple reflections between the sample surfaces, broadens the resonances, 
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and thereby eliminates the spurious effects in the contribution from resonances. In the ex- 
periments, variations in sample thickness break the coherence of multiple reflections and 
eliminate resonance effects. Our use of a damping factor is a simple way to mimic this 
phase breaking and eliminate resonance effects. We use the same damping for all propa- 
gating waves. We also choose this damping to be much weaker than the damping of any 
of the evanescent modes. As a consequence, the damping we introduce reduces the overall 
intensity of the transmitted light without affecting the relative transmission amplitudes of 
different Fourier components of the transmitted light. Thus, the damping that we include 
does not affect the transmission line scans we calculate, other than to change the absolute 
scale of the transmission intensities. We get the same transmission line scans, except for the 
overall scale factor, when we eliminate spurious effects of the resonances by increasing the 
wavevector cutoff. 



To model transmission-NSOM images, we need to calculate the amount of the transmit- 
ted light that is collected by optics with a numerical aperture NA. We use the following 
simple model (see Fig. |l|). We assume that the optics collects all of the transmitted flux 
that passes through an x-y plane in the far- field away from the NCG. All flux leaving the 
NCG within the numerical aperture of the optics, the angle 9 in Fig. [I], is collected. 

To calculate the total intensity of the collected light, J, we first need the time-averaged 
Poynting vector for a harmonic free-space field with frequency u, 



D. Image formation 




(33) 



The flux passing through an x-y plane in the far-field is (S z (x, y, z — > oo)). 



Using a 2D Fourier expansion to write the fields, 



E(x, y,z) = E ( k x, k y ) exp(i(k x x + k y y)) exp(ik z (z - z s )) 



(34) 
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where z s is the position of the back sample- surface, the sum is over all 2D fc-space, k z 
is obtained from the free-space dispersion, and the phase of k z is chosen to give forward- 
propagating or damped waves going into free space away from z s . A similar equation holds 
for B. The integrated flux passing through the x-y plane at z is the total intensity, 



hot = J dxdy (S z (x,y,z)) 



c 

8^ 



Re{ J2 z • (E(^, k y ) x B*(k x , k y )) exp(i(k z - k*)(z - z s ))}. (35) 



For propagating modes k z = k*; for evanescent modes k z = —k*. In the far-field limit, the 
exponential factor vanishes for evanescent modes and is 1 for propagating modes, thus 

J tot = Re{E x (k x , k y )B*(k x , k y ) - E y (k x , k y )B*(k x , k y )}. (36) 

The optics collects light within the cone defined by the numerical aperture, reducing the 
cutoff for light collection from the free-space cutoff 2n / \ to 2irNA/X where NA = sin 6 in 
free space. Thus, the total intensity of transmitted light / collected by the optics is 



c 

8^ 



R e {E x (k x , k y )B*(k x , k y ) — E y (k x , k y )B*(k x , k y )}. (37) 



y/k%+ky< — 



Simulated images are obtained by calculating / as a function of the relative positions of 
the tip and sample. In this paper we present one-dimensional line scans of these images. 



E. Nanochannel glass model 

As shown in Fig. [I], the NCG sample is a 2D triangular lattice of glass rods in a glass 
matrix. We assume that rods are perfect cylinders. In the experiments, the lattice spacing 
is 1.07 /im and the core radius r c = 0.37 fim. We use these values to model the measured 
images. We do additional calculations with other lattice parameters to further test the 
potential of NSOM for probing these structures. The measured samples are approximately 
250 fim thick. We use this thickness to model the measured images. Our results for thick 
samples are insensitive to sample thickness because we eliminate Fabry-Perot oscillations by 
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use of the damping factor. We also consider thin samples where sample evanescent modes, 
excited by the tip field, can be transmitted through the sample. We assume that the index 
of refraction is real and uniform in each glass, has the bulk value in each glass, and has 
a step discontinuity at the core/matrix glass interface. This is the only assumption that 
we make about the dielectric structure of the NCG. The matrix glass has a bulk index of 
refraction of 1.678 for A = 488 nm and 1.657 for A = 670 nm [|7|. The core glass has a higher 
index so it acts as a guiding region. For A = 488 nm, the index difference is An ~ 0.011 
0. For A = 670 nm, the index difference is An ~ 0.019 @. We use these values to model 
the measured images. Other calculations are done with higher An for additional tests. The 
scanned NCG-surface was polished. The core glass etches preferentially during polishing, 
leaving 3.5 ±0.5 nm depressions centered on the channel glasses. The experimental scans are 
done at constant separation, z ap , between the tip aperture and the NCG. For the calculations, 
we model a flat surface, ignoring the effect of these small depressions on the fields at the 
sample surface. Our results show that the transmission NSOM images are determined by 
the coupling of the tip fields to the bulk photonic modes of the NCG and should not depend 
on the surface structure. Also, our results do not change much for nm-scale changes in z ap 
or sample thickness so we can ignore topographic effects on the calculated images. 



F. Significance of NA and fc-space cutoffs 

Our NCG experiments 0, as well as other experiments [18|, show that collection optics 
NA has a strong effect on transmission NSOM image contrast. The collection optics NA 
determines which transmitted modes are collected. Figure |3| shows the reciprocal-lattice 
points for the sample studied experimentally. Also shown are the boundary for the first 
Brillouin zone and the cutoffs for the transverse wavevector that determine which transmit- 
ted modes are collected by the optics for different A and NA (see Eq. [J7|). The transverse 
wavevectors for transmitted modes that are collected by the optics lie inside the cutoff. As 
NA increases, the cutoff for mode components that are collected by the far-field detector 
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increases. Also, for fixed NA the cutoff increases with decreasing A. For NA = 0.28, the 
cutoffs lie well inside (at) the boundary of the first Brillouin zone for A = 670 nm (488 nm). 
For NA = 0.7, the cutoff extends to the edge of the second Brillouin zone for A = 670 nm 
and into the fourth Brillouin zone for A = 488 nm. From our results we find that the scans 
for different A are similar if the scans have similar cutoffs and thus have contributions from 
similar sets of transverse wavevectors. Thus the scans for A = 670 nm are similar to those 
scans for A = 488 nm done at smaller NA. 

IV. RESULTS 
A. Line scans 

First, we show that the theory is able to describe the images of the NCG made with 
transmission NSOM. Calculated line scans along x of the transmitted intensity are shown 
in Figs. f| and |5| for A = 488 nm and 670 nm respectively. The intensities are given in 
arbitrary units. However, the same units are used for these two figures and for all other 
figures that deal with the same NCG, except where it is noted that the results have been 
rescaled. Thus, intensities plotted in different figures for the same NCG can usually be 
compared directly. The dependence on NA is shown. The calculated line scans reproduce 
the key features in the measured line scans. For small NA, the line scans are peaked at 
the channel glass centers, with a peak/valley intensity ratio of about 1.5, even though the 
index contrast is only 1-2%. As NA increases, more of the transmitted light is collected 
and the average transmitted intensity increases. The peaks broaden and flatten and the 
peak/valley ratio decreases. At NA = 0.5 (0.7) for A = 488 (670) nm, the peaks about core 
centers develop structure with depressions in the middle and side peaks ^0.2 /im from the 
centers. These side peaks produce the ring structure observed in the experimental images. 
The experimental images at 0.7 NA show the rings prominently for A = 488 nm but only 
hint that rings are starting to appear for A = 670 nm. Our calculations also show that the 
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rings form at lower NA for shorter A. This occurs because more of fc-space is collected, for 
a given NA, for shorter A. Scans along other directions (see Fig. |6] for scans along y) give 
similar results. In particular, the structure inside the glass core is identical for the scans 
along different directions, confirming the cylindrical symmetry of this structure. This gives 
evidence that the structure in the peak is due to transmitted light that propagates through 
the sample in photonic modes that are concentrated in the cylindrically symmetric cores. 

Figure |2| shows a comparison of the experimental and calculated line scans along x at 
NA = 0.7 for both A. The curves have been scaled to have the same average value and 
then shifted for clarity. There is strong qualitative agreement between the data and the 
calculations, as we have already discussed. The ring structure occurs at similar positions, 
well inside the edge of the channel glass, in both the experimental and calculated scans. This 
ring structure is not directly related to the channel glass edge but is connected to the spatial 
distribution of NCG photonic modes. The most noticeable differences are in the intensity 
ratios. The depressions in the peaks for 488 nm are more pronounced in the experimental 
line scans than in the calculated scans. The experimental data typically shows more contrast 
than the calculations. For example, the small dip in intensity at the three-fold site for a y 
scan is harder to see in the calculated scans (see Fig. |^) than in the experimental data [|7J. 
The experimental intensity ratios at 0.28AM. for peak/three- fold site and bridge/three-fold 
site are 1.58 and 1.18 for A = 488 nm and 1.60 and 1.10 for A = 670 nm |7|. The calculated 
ratios at 0.3NA are 1.69 and 1.05 for A = 488 nm and 1.51 and 1.02 for A = 670 nm. At 
0.7NA the experimental peak/three- fold site intensity ratios are 1.15 for A = 488 nm and 
1.28 for A = 670 nm 0. The calculated ratios at 0.7NA are 1.08 for A = 488 nm and 1.14 
for A = 670 nm. In all cases, the experimental and calculated intensity ratios are similar. In 
all cases except one, the experimental ratio is larger than the calculated ratio. This suggests 
that the calculations may underestimate slightly the coupling to the NCG modes that are 
concentrated in the glass cores and produce the structure in the line scans. Some of the 
small disagreement between experimental and calculated intensity ratios can probably be 
removed by using better models for the NCG and its dielectric profile, the tip field, or the 
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multiple scattering between the tip and the NCG. We have not yet tested these possibilities. 

Consider the following simple analysis of the transmitted intensity to understand the 
structure in the line scans. The index difference between the core and matrix glasses in 
the NCG is small, so much of the collected light will be transmitted through the NCG in 
extended, nearly uniformly distributed modes. The channel glass does have the higher index, 
so some of the collected light will be transmitted in modes more concentrated in the glass 
cores. The line scans can be better understood by separating the collected intensity / into 
these two contributions, I = I m + AI. I m is the minimum intensity. This is the uniform 
intensity that is transmitted and collected, regardless of the tip position. AI is the variation 
in the line scan above the minimum intensity. AI includes the structure in the line scan. 
Figs. f|-| show that I m is approximately the collected intensity at the bridge site, which 
varies approximately as NA 2 . The uniform response of the NCG is proportional to the total 
phase space area of the collected light. The amplitude for variations of AI for a line scan 
along x is the difference in intensity, I c — 1^, between center and bridge sites. As a function 
of increasing NA, I c — lb first rapidly increases, peaks, and then slowly decreases. The 
amplitude for variations in AI is a maximum for line scans collected with NA pa 0.3 — 0.35 
for A = 488 nm and NA pa 0.4 — 0.45 for A = 670 nm. For both A, the maximum amplitude 
for variations in a line scan occurs for collection with X/NA pa 4r c , where r c is the radius 
of the channel glass. The minimum wavelength, X m i n , for the in-plane variation of the 
light collected by the optics for a particular NA is X/NA. The maximum amplitude for 
variations in a line scans occur when the scans are collected with an NA that just accepts 
light concentrated in modes that can fill the channel glass core with half a wavelength (i.e. 
the lowest order "channel" mode with transverse wavelength Xt ran s ~ 4r c pa X m i n , as one 
would expect if the channels were waveguides with strong confinement of the fields). 
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B. Importance of the tip field 



Calculations have been done to determine the effect of structure in the tip field on the 
line scans. In the Bethe-Bouwkamp model for the aperture field, the polarization of E ap is 
determined by the polarization of the plane-wave field incident on the hole in the screen. We 
tried different polarizations of the incident field in the BB model and found no significant 
differences in the calculated line scans. This is consistent with the experimental observation 
that the images were insensitive to polarization. This insensitivity is another indication that 
the local variations in the images arise from the NCG photonic modes that are concentrated 
in the channel-glass cores and have cylindrical symmetry in the cores. 

The effect of the separation, z ap , of the aperture from the NCG and of the aperture 
radius, a, on the line scans has been tested as well. Figure ^ shows the dependence on z ap 
for A = 488 nm and 0.7NA. Similar results were seen for other NAs and for A = 670 nm. As 
the tip/sample separation increases, there is a steady decrease in the total collected intensity 
but little change in the structure in I until the intensity begins to increase and the contrast 
reverses at large z ap [z ap « 500 nm). This suggests that the NCG modes that transmit 
light through the thick film do not couple efficiently to the strongly-localized evanescent 
modes emitted by the tip and that the structure in the line scans is related to the mode 
density in the NCG rather than the field distributions near the tip. Line scans calculated 
by including only the propagating tip modes are almost identical to line scans calculated 
with both the propagating and the evanescent tip modes. This confirms the experimental 
observation that the evanescent tip modes did not couple effectively to the modes in the 
NCG that was imaged. I increases at large z ap because the tip field spreads laterally enough 
to couple efficiently to multiple channel cores at the same time. When z ap is small, the 
tip field couples to one core at a time. The contrast reverses at large z ap possibly because 
the coupling to multiple channels simultaneously is more efficient when the tip is between 
channels rather than centered on a channel. Contrast reversal at large z ap was also observed 
experimentally. 
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The dependence of the line scans on tip radius a is shown in Fig. |^ at 0.7 NA and 
A = 488 nm for the tip close to the NCG, i.e. z ap = 10 nm. Similar results were seen for 
other NA and A. The intensity of the aperture field increases rapidly with increasing a. 
We eliminate this dependence by scaling the line scans to have the same average intensity. 
In the figure, the scans have been shifted for clarity. For small tips with a ~ 100 nm, the 
structure is independent of tip size. Again, this suggests that the NCG modes that transmit 
light through the thick film do not efficiently couple to the strongly localized tip modes. 
The structure persists but weakens for 100 nm ~ a ~ 300 nm and the contrast reverses 
for a > 300 nm. The structure weakens and reverses contrast when the tip becomes wide 
enough to couple to multiple channels simultaneously. 



C. Significance of photonic crystal optical modes 

In this section we show that the structure in the line scans is related more closely to the 
density of the photonic modes in the NCG that are concentrated in the cores and not as 
closely to the density of the more extended, nearly uniform modes or to the structure in the 
tip field. We present line scans for different channel-glass sizes and different NCG lattice 
constants in Figs. § and |10[ The calculations are done for 0.7 NA and A = 488 nm for a 50 



nm tip radius. As the size of the channel glass increases, the width of the central peak in 
the line scan increases. The separation between the side peaks that produce the rings in the 
images also increases. For each core size, the ring structure occurs well inside the channel 
glass and does not give a direct measure of the channel-glass size. Additional structure 
develops as the glass core increases in size and additional modes can be concentrated at the 
cores. The ring structure is eliminated by reducing the core size and making it difficult for 
any mode except the lowest order mode to be concentrated on the cores. The line scans 
show a similar dropoff near the glass-core edge for different core sizes, indicating that the 
spread of core-like modes into the matrix glass, which depends weakly on core-glass size, 
determines the dropoff. The core-glass size and the lattice constant can be varied separately. 
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Fig. ^ shows that the structure in the line scans is determined by the NCG modes that are 
concentrated to the glass cores. Increasing the NCG lattice constant increases the spacing 
between the main features in the scan without changing the structure in those features. 



D. High-contrast samples 

The structure in the experimental images and in the line scans discussed so far results 
because the NCG supports photonic modes which are concentrated on the channel-glass 
cores. The structure is insensitive to the tip field when the tip field only couples to one 
channel at a time. Contrast reversal occurs when the tip field couples to multiple channels 
simultaneously. Transmission of the evanescent near fields from the tip through the NCG 
is negligible in the samples studied for two reasons. First, the transverse wavevector k ev 
for localized evanescent modes emitted by the tip is much larger than the reciprocal-lattice 
basis vectors of the NCG samples, i.e. k ev « it /(2a) ^> G®,G®. The evanescent modes can 
couple to the NCG photonic modes in the experimental samples only through very high- 
order Fourier components of the index-of-refraction profile. To enhance this coupling, one 
can reduce the NCG lattice constant so that the tip field and the index profile vary on the 
same length scale in the plane, i.e. so that k ev is comparable to G° x and G®. One can also 



increase the index contrast An. In Fig. [11] we show line scans calculated for an NCG with 
half the lattice constant and half the core-glass radius of the experimental samples, with 
the same index for the matrix glass as in the experimental samples, but an index contrast 
An = 1 which is much bigger than in the experimental samples. The sample thickness 
was chosen to be 250 /iin as in the experimental samples. The calculations were done as a 
function of NA for A = 488 nm with a 50 nm radius tip, scanning 10 nm from the NCG. The 
structure in the line scans is similar to that obtained for the experimental samples with only 
1-2 per cent index contrast. At low NA, a single central peak occurs. As NA increases, this 
peak broadens, flattens and then evolves into the ring structure. However, the peak/bridge 
site contrast ratios are much larger when An = 1 . A comparison of line scans of the high An 
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sample done at z ap = 10 nm with line scans done at z ap = 50 nm (see Fig. |T^) shows a much 
stronger dependence on z ap than seen for low An samples (Fig. |7|). There is a much larger 
decrease in total transmitted intensity for the same change in tip height when An = 1. More 
importantly, the structure in the line scans of the high An sample changes significantly as 
the tip height changes from 10 nm to 50 nm. The ring structure disappears when z ap is 
increased from 10 nm to 50 nm above the high An sample. 

These calculations were done for a 250 /zm thick sample. Many of the NCG modes that 
couple to tip evanescent fields are evanescent modes. These NCG modes make a negligible 
contribution to transmission through a thick sample. Transmission of evanescent fields from 



the tip through the sample can be enhanced by reducing the sample thickness. Figure 13 



shows line scans done for the same sample as for Fig. |TT] except that the sample thickness 
is only 0.1 /im. Increased transmission of the evanescent tip fields by the evanescent sample 
modes greatly enhances the structure in the line scans. The ring structure is much more 
prominent. Additional structures at interstitial sites, such as the bridge sites, is now present. 
A comparison (Fig. ITJ) of the line scans for z ap = 10 nm and z ap = 50 again shows a strong z ap 
dependence both at high and low NA. Calculations for samples which are approximately 
a wavelength thick show additional Fabry Perot structure. The thickness of the sample 
determines which sample modes provide enhanced transmission. The tip/sample separation 
determines how strongly these modes are excited. 



V. CONCLUSIONS 

A theory has been presented to model recent transmission NSOM images made of 
nano channel-glass arrays. The theory describes the entire process of image formation, in- 
cluding the tip field that is the light source, the photonic modes of the nano channel-glass 
array, coupling of the tip field to the sample modes, and transmission of the fields to the 
collection optics. The theory is able to reproduce the key features observed in the experi- 
mental images, including the evolution of the images with increasing numerical aperture of 
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the collection optics and the dependence on wavelength. Most importantly, the theory is 
able to explain a ring-like structure that develops in the images as the numerical aperture 
is increased. 

For the samples studied experimentally, transmission NSOM probes most directly the 
photonic mode structure of the nano channel-glass array. The structure in the images can 
be related to the photonic modes of the array that are strongly concentrated at the glass 
cores. The numerical aperture of the collection optics can be varied to control the size of the 
phase space of the photonic modes that is sampled by the transmission experiments. Angle- 
resolved mapping of the photonic modes could be achieved by collecting light transmitted 
at specific angles. The structure in the images is insensitive to the details of the tip field 
because the samples do not couple strongly to the localized evanescent modes emitted by 
the tip. In these experiments the images become sensitive to the tip field only if the tip field 
is so broad that it can simultaneously excite multiple channel-glass cores in the array. 

Contrast and structure in the images can be enhanced greatly by increasing the transmis- 
sion of the evanescent tip fields through the sample. This can be accomplished by choosing 
samples with higher index contrast or with index profiles with spatial variations that more 
closely match the spatial localization of the tip field. This enhanced contrast comes at the 
price of increased sensitivity to the tip field. This increased sensitivity can be exploited to 
better probe the sample, but it must be understood to properly separate sample properties 
from probe properties. 
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FIGURES 

FIG. 1. Schematic of the transmission NSOM experiment. The nanochannel-glass array sample 
is scanned at constant separation across the NSOM tip. The transmitted light is collected by an 
objective with numerical aperture NA. Key elements of the theory are represented schematically. 
The tip field is described by the Bethe-Bouwkamp model |IO| , pj]1 . The collected light is modeled as 
the total flux that leaves the sample inside the acceptance cone defined by AM and passes through 
a plane, as indicated, in the far field. The array geometry is indicated. 

FIG. 2. Transmission as the tip scans along x: experimental scans (□,+) and calculated scans 
( — ) scaled to have the same average intensity and shifted for clarity. The core-glass centers are 
marked by vertical dashed lines. At x = the tip and core glass centers are aligned. 

FIG. 3. Reciprocal lattice points (•) and first Brillouin zone boundary ) for the NCG 
sample. Cutoff wavevectors for A = 488 nm ( — ) and 670 nm (• • ••) are shown. Each inner (outer) 
circle is for NA = 0.28 (0.7). 

FIG. 4. Transmitted intensity line scan along x: A = 488 nm, tip/sample separation 
Zap = 10 nm, tip radius a = 50 nm. Dependence on collection optics NA is shown. Positions 
of channel-glass centers (c) and bridge sites (b), the midpoints on the lines joining adjacent glass 
centers, are shown. 

FIG. 5. Transmitted intensity line scan along x: A = 670 nm, z ap = 10 nm, a = 50 nm. 
Dependence on NA is shown. Positions of channel-glass centers (c) and bridge sites (b) are shown. 

FIG. 6. Comparison of transmitted intensity line scans along x ( — ) and y (• • ••): A = 488 nm, 
z ap = 10 nm, a = 50 nm. Positions of channel-glass centers (c), bridge sites (b) along x, and 3-fold 
sites (3f)along y are shown. The 3-fold symmetric sites are at the middle of the triangles formed 
by 3 adjacent nearest-neighbor glass centers. 

FIG. 7. Transmitted intensity line scan along x: A = 488 nm, a = 50 nm, NA = 0.7. Depen- 
dence on tip/sample separation, z ap , is shown. 
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FIG. 8. Transmitted intensity line scan along x: A = 488 nm, z ap = 10 nm, NA = 0.7. 
Dependence on tip radius, a, is shown. The line scans have been scaled to have the same average 
intensity (1000 on the scale used) and then shifted for clarity from the scan for a = 100 nm by 
multiples of 100. 

FIG. 9. Transmitted intensity line scan along x: A = 488 nm, z ap = 10 nm, NA = 0.7, 
a = 50 nm. Dependence on channel-glass core radius is shown. 

FIG. 10. Transmitted intensity line scan along x: A = 488 nm, z ap = 10 nm, NA = 0.7, 
a = 50 nm. Dependence on lattice constant and channel-glass core radius is shown. The curves 
are all on the same scale but have been shifted for clarity by from the scan for the experimental 
sample (lattice constant 1.07 /im and core radius r c = 0.37 um) by multiples of 300. The positions 
of the core centers (c) and bridge sites (b) for the NCG studied experimentally are indicated. 

FIG. 11. Transmitted intensity line scan along x: A = 488 nm, z ap = 10 nm, a = 50 nm. 
Dependence on NA is shown. The NCG lattice constant is 0.535 ^m, the channel-glass core radius 
is 0.186 /im, and the index contrast is An = 1.0. The NCG sample is 250 /im thick. The vertical 
lines indicate channel-glass centers (c), bridge sites (b), and edges of the channel glass (e). 

FIG. 12. Comparison of transmitted intensity line scan along x: A = 488 nm, a = 50 nm, 
z ap = 10 nm ( — ) and 50 nm (• • ••). The NCG has the same lattice parameters, thickness, and 



index contrast as in Fig. 11. Dependence on NA is shown. 



FIG. 13. Transmitted intensity line scan along x: A = 488 nm, z ap = 10 nm, a = 50 nm. 
Dependence on NA is shown. The NCG lattice constant is 0.535 fim, the channel-glass core radius 
is 0.186 fim, and the index contrast is An = 1.0. The NCG sample is 0.10 fim thick. 

FIG. 14. Comparison of transmitted intensity line scan along x: A = 488 nm, a = 50 nm, 
Zap = 10 nm ( — ) and 50 nm (• • ••). The NCG has the same lattice parameters, thickness, and 



index contrast as in Fig. 13. Dependence on iVj4 is shown. 
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